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Abstract 

We present a new approach to constructing inflationary models in 
closed universes. Conformal embedding of closed-universe models in a 
de Sitter background suggests a quantisation condition on the available 
conformal time. This condition implies that the universe is closed at no 
greater than the 10% level. When a massive scalar field is introduced to 
drive an inflationary phase this figure is reduced to closure at nearer the 
1% level. In order to enforce the constraint on the available conformal time 
we need to consider conditions in the universe before the onset of inflation. 
A formal series around the initial singularity is constructed, which rests 
on a pair of dimensionless, scale-invariant parameters. For physically- 
acceptable models we find that both parameters are of order unity, so 
no fine tuning is required, except in the mass of the scalar field. For 
typical values of the input parameters we predict the observed values of 
the cosmological parameters, including the magnitude of the cosmological 
constant. The model produces a very good fit to the most recent CMB 
(cosmic microwave background) radiation data. The primordial curvature 
spectrum provides a possible explanation for the low-£ fall-off in the CMB 
power spectrum observed by WMAP (Wilkinson anisotropy probe). The 
spectrum also predicts a fall-off in the matter spectrum at high k, relative 
to a power law. A further prediction of our model is a large tensor mode 
component, with r ~ 0.2. 

PACS numbers: 98.80.Cq, 98.80.Jk, 98.80.Es, 04.20.Gz 

1 Introduction 

Experimental evidence for the existence of a cosmological constant has grown 
dramatically over recent years. It now seems likely that the cosmological con- 
stant, or some form of 'dark energy', is responsible for around 70% of the total 
energy density of the universe . There are two attitudes one can take towards 
the dark energy. One is that it is a relic from symmetry breaking in a more 
complete theory, and as such has a field-theoretic origin. As is well known, 
this idea runs into difficulties in explaining the current magnitude of the cos- 
mological constant. The second viewpoint is the more geometric one, that the 
cosmological constant is related to the large scale geometry of the universe. This 
is favoured by the fact that the energy density in the cosmological constant is 
currently of the same order as that of the large scale matter density. But this 
viewpoint is silent on any mechanism that might explain the size of A. Here we 
develop the geometric viewpoint, and suggest a novel boundary condition that 
can account for the size of the cosmological constant. 
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A two-dimensional space of constant negative curvature and positive-definite 
signature can be conveniently visualised using the Poincare disk construction |2 
[3]. The disk represents the entire 2-dimensional space, with geodesies rep- 
resented by circle arcs that intersect the boundary at right angles. Here we 
generalise this idea to provide a Lorentzian picture of de Sitter space. Open, 
closed and flat cosmological models can be embedded into sections of this space. 
The case of closed models is of greatest interest here, and for these we see that a 
natural boundary condition presents intself. This places the initial singularity at 
the midpoint of the conformal picture of de Sitter space, which implies that the 
total amount of conformal time available to the universe should equal tt/2. For 
a simple dust cosmology this boundary condition picks out a single trajectory in 
the £Im~Qa plane, and for the present epoch predicts a universe that is closed 
at around the 10% level. When coupled with knowledge of £Im, the boundary 
condition therefore predicts a value of £Ia that is close to the experimental value, 
in sharp contrast to predictions of A from high energy physics. 

To further reduce the predicted value of Sltot we next consider an inflationary 
phase. This will still be necessary to seed the initial perturbations, and has 
the added benefit of using up some conformal time before the universe enters 
its matter-dominated phase. We restrict attention to the simplest case of a 
single, massive scalar field <f>, which ensures that the field equation for 4> is 
linear. This is the most economic model that involves the least new physics. In 
order to apply our boundary condition we need to study the evolution of the 
scalar field from the initial singularity, through the inflationary region, before 
matching onto a standard cosmological model. Expanding the field equations 
around the initial singularity involves a complicated iterative scheme that can be 
extended to arbitrary precision. The series is governed by two parameters, which 
effectively control the degree of inflation and the curvature. Fixing both of these 
parameters to be of order unity produces inflationary models in a closed universe 
that are consistent with observation. This runs counter to the claim that it is 
difficult to obtain closed universe inflation without excessive fine tuning |SJ 

ElCllHj. 

An attractive feature of closed-universe models is that one always has a fixed 
distance scale to refer to, fixed by the size of the universe at the time of interest. 
One can use this, for example, to fix the onset of the regime in which quantum 
gravity effects might be expected to dominate. This will be when the radius of 
the universe is of the order of the Planck scale. For the regions of parameter 
space of interest here, we find that quantum gravity effects are not relevant until 
well before the onset of the inflationary region. We therefore have little choice 
but to run our equations back in time to close to the initial singularity. This 
approach has the advantage of setting the values of the fields as the universe 
enters the infaltionary regime, making the model highly predictive. 

As the universe exits the inflationary region it evolves as if it had started from 
an effective big-bang, with a displaced time coordinate. By this point a sizeable 
fraction of the available conformal time has been used up, which further restricts 
the extent to which the universe can deviate from flatness. The result of these 
effects is the imposition of a see-saw mechanism linking the current state of the 
universe and the initial conditions. The more we increase the number of e-folds 
during inflation, the smaller the value of the cosmological constant, and vice- 
versa. Furthermore, we find that A scales as exp(— 6N) 7 where N is the number 
of e-foldings during inflation, and is in the range of 40 - 50. The smallness of A 
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is then an immediate consequence of this model. With initial conditions chosen 
to give the required number of e-foldings to generate the observed perturbation 
spectrum, we find that the predicted universe is closed at the level of a few 
percent. Given a present value of Qm at around 0.3 we single out models with 
Oa ~ 0.7. That is, the conformal time constraint predicts both the degreee of 
flatness of the universe and the magnitude of the cosmological constant. 

For this model to be valid it has to correctly predict the growth of pertur- 
bations during the inflationary phase. The primordial curvature spectrum we 
predict shows a strong departure from a simple power law, with an exponential 
cut-off at low wave number, and an additional fall-off at high k. The resulting 
CMB power spectrum is in good agreement with the WMAP data jT], and is 
able to explain the observed low-£ deficit. 

We write the FRW equations as 



R 2 + k A 8nG 

where k = 0, ±1 and throughout we work in units where c = h = 1. The 
scale parameter R has dimensions of distance and for closed models R is the 
proper radius of the universe. This is why we prefer the symbol R over the more 
common a. Where they add clarity, factors of G are included, so that G has 
dimensions of (distance) 2 . The ratios flu and J7a are defined by 

_^G P A 
where H = R/R. The universe is spatially flat if Um + = 1- 



2 Picturing de Sitter space 

A de Sitter space is a space of constant negative curvature. As such, it forms the 
Lorentzian analogue of the non-Euclidean geometry discovered by Lobachevskii 
and Bolyai Ej- Two-dimensional non-Euclidean geometry has an elegant 
construction in terms of the Poincare disc, in which geodesies arc represented as 
'd-lines' — circles that intersect the disc boundary at right-angles (see figure^. 
Here we provide a similar picture for 2-dimensional de Sitter space, which is 
sufficient to capture all of the key features of the geometry [Hj. We start with 
an embedding picture, representing de Sitter space as the 2-surface defined by 

T 2 - X 2 - Y 2 = -a 2 , (4) 

where (T,X,Y) denote coordinates in a space of signature (f,2) and a is a 
constant. The resulting surface is illustrated in figure [21 The figure illustrates 
some key features of the global de Sitter geometry. In particular, spatial sections 
are closed, whereas the timelike direction is open. So de Sitter space describes a 
closed universe that lasts for infinite time. One can also set up local coordinate 
patches for which spatial sections are flat or open, but these coordinates are 
not global and only cover part of the full manifold. Null geodesies are straight 
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Figure 1: The Poincare disc. Points inside the disc represent points in a 2- 
dimensional non-Euclidean (hyperbolic) space. A set of (i-lines are also shown. 
These are (Euclidean) circles that intersect the unit circle at right angles. Given 
a d-line and point A not on the line, one can find an infinite number of lines 
through A that do not intersect the line. 



lines formed by the intersection of the surface and a vertical plane a distance a 
from the timelike axis. Despite the fact that the space is spatially closed, the 
furthest a photon can travel is half of the way round the universe. 

To establish an analogous picture to the Poincare disc for de Sitter space, 
we start by considering the spatial section at T — 0. This section is a ring of 
radius a, which is mapped onto a straight line in a (1, 1) Lorentzian space via a 
stereographic projection. Null geodesies from this section are now represented 
as 45° straight lines in Lorentzian space. Since geodesies from opposite points 
on the ring meet at infinity, we arrive at a boundary in the timelike direction 
defined by a hyperbola. This construction provides us with a Lorentzian view 
of de Sitter geometry. Timelike geodesies in de Sitter space are represented 
by hyperbolae that intersect the boundary at a right-angle (see figure EJ). Here 
'right-angle' is defined in its Lorentzian sense (the tangent vectors have vanishing 
Lorentzian inner product). This view has the convenient feature that geodesies 
are lines of constant Lorentzian distance from some fixed point. The metric 
associated with this view of de Sitter geometry has the conformal structure 

ds " = (fl 2 +x a 2 4 _ t2)2 ( rft2 - dx2 )> (5) 

which makes it clear that null geodesies must remain as straight lines in the x-t 
plane. 

There are many fascinating geometric structures associated with de Sitter 
geometry, many of which mirror those of non-Euclidean geometry. For exam- 
ple, one can always find a reflection that takes any point to the origin jlllj . One 
can then prove a number of results at the origin, where the geodesies are all 
straight lines, and the results are guaranteed to hold at all points. A similar ap- 
proach can be applied to anti-de Sitter space, except now the diagram is rotated 
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Figure 2: Two-dimensional de Sitter Space. The timelike direction is vertical, 
and spatial sections are closed. The right-hand diagram shows a null geodesic, 
which is a straight line in the embedding space. 




Figure 3: The Lorentzian view of two-dimensional de Sitter space. The bound- 
ary is defined by two hyperbolae (shown with thick lines). All geodesies through 
the origin are straight lines, and null geodesies are always straight lines at 45°. 
Two further geodesies, one spacelike and one timelike, are also shown. These 
are hyperbolae which do not pass through the origin. The timclikc geodesic in- 
tersects the boundary in such a way that the two tangent vectors have vanishing 
Lorentzian inner product. 
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Figure 4: Flat and open cosmological models, de Sitter space contains sections 
representing both flat and open spacetimes. A flat spacetime (left) consists 
of the space contained inside a light-cone located on past infinity. An open 
spacetime fits inside the lightcone from the origin. Both pictures can be used 
to illustrate cosmological models, with the initial singularity represented by a 
hyperbolic spacelike surface. An example of these surfaces is shown on each 
diagram. 

through 90°, as it is the timelike direction that is formed from a stereographic 
unwrapping of a circle. 

3 Cosmological models and the de Sitter end 
state 

All cosmological models are conformally flat, and can all be interchanged via 
conformal transformations. Furthermore, all models containing a cosmological 
constant, and which do not recollapse, will tend towards a de Sitter end. Such 
models should fit neatly within the conformal diagram of figure [3] with future 
infinity represented by the upper hyperbolic boundary. There are three cases to 
consider, for open, closed and flat cosmologies. In this paper we are interested in 
closed cosmologies, but before looking at these it is informative to consider the 
other cases. A flat section of de Sitter space corresponds to a region contained 
within a light cone from a point located on the boundary at past infinity (see 
figure UJ. Surfaces of constant cosmic time are represented as hyperbolae, any 
one of which can then be chosen to represent the initial singularity. Similarly, 
an open section of de Sitter space is represented by the area inside a light-cone 
from a point in the middle of the de Sitter picture. A more general open A- 
cosmology will have the initial singularity located on a spacelike hyperbola. (See 
Ratra and Peebles m for more details of inflation in open cosmologies.) 

The diagrams for flat and open cosmologies make it clear that one is not 
employing the full de Sitter geometry in a symmetric manner. Closed models, 
on the other hand, have a more natural embedding, with the initial singularity 
represented by a spacelike surface in the conformal diagram of figure 01 Here we 
begin to see a natural boundary condition emerging. Suppose that we insist that 
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t = represents the initial singularity, so that the entire history of the universe 
can be conformally mapped into the top half of the de Sitter diagram. (This does 
not imply that the actual universe is static, it is merely a statement about its 
conformal representation.) This is certainly the most natural, symmetric place 
to locate the singularity. An immediate consequence is that any photon emitted 
from the initial singularity may travel a maximum of precisely one-quarter of 
the way round the universe in the entire future evolution of the universe. An 
alternative way of saying this is that the past horizon projected back to t = 
should cover half of the de Sitter geometry. 

Given an equatorial angle on a 3-sphere, a photon travelling round the 
equator will satisfy 

d(f> 1 , s 

* = fl' (6) 
where t is cosmic time. Traversing the entire universe corresponds to running 
from < 4> < 2n. If a photon is only to travel one-quarter of the way round we 
therefore require that 

00 dt _ f 00 dR 

(l R~J R(-l + AR 2 /3 + 87^/3) V2 ~ 7r/ ' 1,1 

At this point is is useful to define the conformal time r\ by 

f* dt ' 

Note that for closed models the conformal time r\ is not the same as the time- 
like conformal coordinate t (see appendix |A"|> . Our boundary condition now 
states that the end of the universe corresponds to a total passing of tt/2 units 
of conformal time. 

The constraint 77(00) = tt/2 can be arrived at through an alternative route 
that works entirely within the conformal representation of cosmological models, 
with the line element written in the form 

ds 2 = -p (dt 2 - dx 2 - dy 2 -dz 2 ). (9) 

The precise form of / depends on the spatial topology. In appendix we show 
that the three cases reduce to 

f = , \ = Q I -rs — n\ t" closed 

J Rsin[rj) J \\ 2 + r 2 -t 2 J A 

o \ 

/ = — flat (10) 

t 1 ~( 2At \ t 

f= Rl^Mv)= 9 { \ 2 +t 2 -r 2 ) A ° Pen 

where r 2 = x 2 + y 2 + z 2 , A is a constant with dimensions of distance, and the 
forms of g and g depends on the model for the matter content. The coordinate 
ranges are as illustrated in figures|21and0J A de Sitter space centred on t = = t 
has the conformal line element 

1 r>\ 2 

ds 2 = K{x2 l l X r2 _ t2) (dt 2 dx 2 - dy 2 - dz 2 ) . (11) 
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Clearly the only solutions that have any chance of matching onto this final state 
are those for a closed universe. Furthermore, the function g must satisfy 



/ \2 A \ 1/2 i 

lim«Kx)= — - (12) 

X^°o V 3 / Y 



where 



2At 



- A2 

But since 



= tan(r?). (13) 



9 = R^f) (14) 

we must then have Rcos{rf) tending to a constant at large times. This is only 
possible if r\ tends to 7r/2, recovering our earlier boundary condition. This 
derivation is instructive in that it reveals how the constraint can be imposed as 
a straightforward boundary condition on a differential equation. In this case, 
the equation for g(x) is 

2,, , 2^f d 9V, d . 2 8irG\ 2 p A 2 A 

The task then is to solve this equation subject to the boundary condition that 
g falls off as l/x for large \- Viewed this way the constraint can be thought of 
as a 'quantisation condition' applied as the universe is formed, which one might 
hope would emerge from a quantum theory of gravity. 

In order to understand the implications of our boundary condition, we first 
write the FRW equations in the form 

dfl A = Qa ((l + 3 7 )n M -2(Q A -1)) 
(«M flM((l + 37)(OM-l)-2n A )' 

where 7 encodes the equation of state in the form 

P = jp. (17) 

Equation 1)16(1 defines a series of flow lines in the 17m _ ^a plane, and also high- 
lights the significance of the points (1, 0) and (0, 1), which are attractors of the 
flow lines. This behaviour is illustrated in figure which shows a series of flow 
lines starting from Qm = 1, = 0. 

An alternative version of the equations, more useful in describing the uni- 
verse since recombination, is to assume that the matter density is made up 
of decoupled dust and radiation. Writing the densities of these as p m and p r 
respectively, we introduce the quantities 

8irGp m 8nGp r 
= 3H 2 ' = 3H 2 ' ' ' 

Since the matter and radiation are decoupled, both stress-energy tensors satisfy 
separate conservation laws, which reduce to 

p m R 3 = constant, and p r R 4 — constant. (19) 
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Figure 5: Evolution curves in the flM~^A plane. The left-hand plot is for dust, 
and the right-hand plot is for radiation. In both cases the curves converge to 
f?A = 1, representing a late-time de Sitter phase. 



For this case we find that the equations governing Qa, and fl r can be solved 
exactly, with the solution governed by two arbitrary constants. We label these 
a and /3, where 

tl 2 m n A (8irG Pm R 3 ) 2 A 

a (n m + n r + n A - 1) 3 27 1 ' 

and 

n _ O^A _ 8irG Pr R 4 A 

ap (n m + n r + n A - 1) 2 9 1 ' 

Present knowledge of f2 m , il r and Q\ fixes a and /3, and so picks out a unique 
curve. The curve will be valid back to last scattering, beyond which the equation 
of state becomes more complicated. 

The diagrams in figure [5] highlight the fact that as the universe evolves the 
flow lines refocus around the spatially flat case, £Im + = 1- One can show 
that, for a large range of initial conditions, by the time we reach the current 
value of fijif s» 0.3 most models are not far off spatial flatness. This predic- 
tion contrasts with models without a cosmological constant, where any slight 
deviation from the critical density in the early universe is scaled enormously by 
the time we reach the present epoch, implying that the parameters in the early 
universe are highly fine-tuned. One can see the problem straightforwardly by 
writing 

K = fl m + O r + n A - 1. (22) 
The parameter n measures the deviation from spatial flatness, and satisfies 

- = #(n m +2fi r -2fi A ). (23) 

K 

If fl A = the right-hand side is always positive, and nothing prevents the 
continued growth of k. The presence of a A term changes this completely. At 
some finite time 20a will overtake the matter terms, and k is refocused back 
to k = 0. The presence of a cosmological constant therefore goes someway to 
solving the flatness problem on its own, without having to invoke inflation. 
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Figure 6: Critical paths as predicted by the de Sitter embedding. The solid line 
represents a matter-dominated universe, and the broken line shows radiation, 
for comparison. The straight line is the critical case of spatial flatness. By the 
time we reach Q,m ~ 0.3 the universe has been driven close to critical density. 

For the case of a universe filled with non-interacting dust and radiation we 
can write the boundary condition of equation Q as 



This equation is therefore an eigenvalue problem for the dimensionless ratios a 
and p. For the case of dust (f3 = 0) the solution turns out to be a ~ 40.468. 
Pinning down a value of a and p picks out a unique trajectory in the Qm—^a 
plane. For the straightforward cases of dust ((3 — 0) and radiation (a = 0) 
we arrive at the two curves shown in figure [5] As the universe is expected 
to be matter dominated for most of its history, the solid line in figure [S] is 
the more physically relevant one. Taking the present-day energy density to be 
around Qm — 0.3, we see that 51a = is predicted to be around J7a ~ 0.83. 
That is, a universe that is closed at around the 10% ratio. Such a prediction 
is reasonably close to the observed value, though it is ruled out by the most 
recent experiments PQ. In order to improve the prediction, we need to use up 
a greater fraction of the conformal time before we enter the matter-dominated 
phase. Such a process is also necessary to solve the horizon problem, and the 
simplest means of achieving this is via an inflationary phase. 

4 Scalar fields and inflation 

Much effort has gone into constructing quantum field theories in de Sitter 
space \M\ I15| - There are two reasons why this is more complicated 

that in Minkowski space: the closed spatial sections give rise to a quantization 
condition; and the presence of event horizons give rise to ambiguities in the 
definition of the vacuum. Similar issues can arise when considering an evolving 




(24) 
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closed universe and will be relevant for a complete treatment of the growth of 
perturbations seeded by a scalar field. In this section we will ignore such issues, 
delaying a discussion of the quantization procedure until section [5] Our model 
for the matter therefore consists simply of a real, time-dependent, homogeneous 
massive scalar field. For a range of initial conditions this system shows the 
expected inflationary behaviour. The equations for this model are 

H + H 2 -\ + ^(2$ 2 ~m 2 tf) = Q (25) 

and 

4> + 3H<j> + m 2 (t) = 0. (26) 

Here <j> has dimensions of (length) , which is the convention widely adopted in 
the cosmology literature (though this does differ from the relativistic quantum 
mechanics literature As stated in the introduction, it is convenient to 

keep in explicit factors of G, to keep the units clear. For closed models the scale 
factor is given explicitly by 

^ = ^ 2 + -V)-H 2 + ^. (27) 

Given a solution to the pair of equations (|25(l and 1(26(1 , a new solution set 
is generated by scaling with a constant a and defining 

H'(t) = aH(at), cf>'(i) = <j>{at), m! = am, A' = a 2 A. (28) 

This scaling property is valuable for numerical work, as a range of situations 
can be covered by a single numerical integration. Furthermore, many physically 
interesting quantities turn out to be invariant under changes in scale. These 
include the conformal time 77, which is easily confirmed to be scale invariant 
from equation (jHJ. Invariant quantities such as r\ turn out to be extremely 
helpful when considering the initial conditions as the universe emerges from the 
big bang. The scaling property does not survive quantisation, however, so has 
to be employed carefully when considering vacuum fluctuations. 

The initial conditions for equations 1(25(1 and ((26(1 are usually set at the start 
of the inflationary period, where they are viewed as arising from some form of 
quantum gravity interaction. But in order to apply our boundary condition we 
need to track the equations right back to the initial singularity, as this is the 
only way that we can keep track of the total conformal time that has elapsed. 
As we argue below, evolving the inflationary equation backwards in time is 
entirely justified, as we do not expect quantum gravity to play a role until much 
earlier in the history of the universe. Inflation on its own does not eliminate 
singularities from cosmology ^7]. Our approach has the added benefit that the 
conditions at the start of inflation are fully determined by a pair of parameters, 
making the model highly predictive. 

As the universe emerges from the big bang the dominant behaviour of H 
is to go as l/(3t). Equation 125(1 then implies that <f> must contain a term 
going as ln(i). But this in turn implies that H must also contain a term in 
iln(i), in order to satisfy equation ((26(1 . Working in this manner we conclude 
that a series expansion in terms of ln(i) is required to describe the behaviour 
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around the singularity. At this point it is convenient to define the dimcnsionless 
variables 

u= v " = ^ (29) 

with t p and m p the Planck time and mass respectively. The series expansion 
about the singularity at t = can now be written 



1 oo 

%) = -^ff,(«)k'W (30) 

^W^-^^Wln'W, (31) 

which ensures that the expansion coefficients are all dimensionless. Substituting 
these into the two evolution equations, and setting each coefficient of ln(u) to 
zero, we establish that 

2 

dHo tt2 uA 8nu ( dfa \ 16ir<fii dfa 



1 du ' 3 3 \ du J 3 du 

jMj , 4tT^ 2 U0 2 

~lu~ + 3 ' [ ' 

with further algebraic equations holding for H2,fa, H^,fa, and so on. So, by 
specifying H , fa and fa, all the remaining terms in the solution are fixed. 
The aim then is to choose the three input functions to ensure that successive 
terms in the series get progressively smaller. This turns out to provide just 
the right number of equations to specify all coefficients, save for two arbitrary 
coefficients in fa. The end result is a series controlled entirely by a pair of 
arbitrary constants: the expected number of degrees of freedom once we have 
fixed the singularity to t = 0. In order to generate curvature it turns out that 
the input functions need to be power series in u 1 / 3 , which ensures that the scale 
factor goes as w 1 / 3 at small times. The series solution is only required to find 
suitable initial conditions for numerical evolution, so only the first few terms 
are required. Expanding up to order u 5 / 3 we find that 

1 32V37T, -,/, 



'2i? A 4tt 2l2 4^3^ 2, \ 66567rfo? 

£ + 3 + Y »X + —»% j u - V 5/ ' ^ 



+ 



to = bo + m 4 ' 8 - H^!^/ 3 



- — I— ( llV3^ 2 - 54V3^A - 216V3ir 3/2 n 2 bl + 36irn 2 b ) u 2 (34) 

12907T V / 



and 

1 M 2 



\/3tt + 36tt6 ) u 2 . (35) 



12tt 216tt 

Under scaling the three key parameters in the model transform as 

6 i-» 6 , bi ^ a 4/3 b4,. (36) 
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Figure 7: Hubble function entry into the inflationary regime. The Hubble func- 
tion emerges from the big bang going as l/3t. As H falls it quickly enters the 
inflationary region, during which it falls linearly. Time t is measured in units 
of t p , and H in units of t~ x . The input parameters were set to bo = 2.48, 
&4 = —0.51 and m — m p (fi = 1). These values of bo and &4//U 4 / 3 are illustrative 
of physically realistic models. The cosmological constant term has no effect in 
this region and can be set to zero. 



The scaling transformation for 64 is entirely as expected, given that it is the 
coefficient of m 4 / 3 in the series for </>q. It follows that the quantity 64 //i 4 / 3 is 
scale invariant. 

The plot in figure illustrates the general behaviour around the initial sin- 
gularity. As the universe emerges from the big bang the energy density in the 
scalar field is dominated by the <p term, and the field behaves as if it is massless. 
It follows that H initially falls as l/(3i). But once H has fallen sufficiently far 
we enter a region in which m 2 (j) 2 starts to dominate over <j) 2 . These are suitable 
initial conditions for the universe to enter an inflationary phase. By varying 60 j 
64 and fi we control both the values of the fields as we enter the inflationary 
period, and how long the inflationary period lasts. The cosmological constant 
plays no significant role in this part of the evolution. The dynamics displayed 
in this figure is quite robust over a range of input parameters. 

The fact that 64 controls the curvature can be seen from equation l|27|) which, 
to leading order, yields 

* 1 ^V'Y-^Y" 0-)-" +.... (37, 



l p [i V 12544717 

Clearly, the arbitrary constant 64 must be negative for a closed universe. The 
terms on the right-hand side of equation 137f) are all scale invariant, apart from 
the first factor of From figure we can see that, for the typical values of 
60 and 64//J, 4 / 3 of interest, the onset of inlation corresponds to \m ~ 0.1. As we 
discuss below, in order to generate perturbations consistent with observation the 
scalar field must have a low mass, with /1 of the order of 10 _6 m p . It follows that 
the onset of inflation occurs at a time of around 10 5 Planck times. Furthermore, 



13 



the radius of the universe at the onset of inflation is given approximately by 

0.2 

R ~ — l p , (38) 
A 4 

and is of order 10 5 Planck lengths. Inflation therefore starts at an epoch well into 
the classical regime. Quantum gravity effects would be expected to be relevant 
when the radius of the universe is of the order of the Planck scale, which occurs 
when u ~ fi 2 and is well before any inflationary period (for physical values of 
/j,). We are therefore quite justified in running the evolution equations back past 
the inflationary regime, and right up to near the initial singularity. It is only 
when R = l p that the equations will break down, and we would look to quantum 
gravity to explain the formation of the initial, Planck-scale sized universe. We 
might hope here that quantum gravity will provide a probability distribution 
for the two classical parameters controlling the evolution of the scalar field. 
Furthermore, we can argue more strongly that we have to run the equations 
backwards in time to well before the start of inflation before reaching an epoch 
where new physics could be expected to enter the problem. 

As the universe is described by a 3-sphere of radius R, the total volume of 
the universe is given to leading order by 

Following from this, an interesting calculation we can perform in a closed uni- 
verse is to find the total energy contained within it in the scalar field. By 
integrating the energy density we find that 

7T / 2187 \ 3/4 1 h /_ M 4/3x3/ 2 n 

We return to a discussion of this equation in sectional 

The typical behaviour as we exit the inflationary region is shown in figure |H1 
The end of inflation is characterised by H w /i. Beyond this point the scalar 
field enters an oscillatory phase, with the time-averaged fields satisfying the 
conditions for a simple matter-dominated cosmology. In this case H can be 
approximated by a curve going as 2/3(i — to), describing a dust model with a 
displaced origin. The universe then appears as if it has been generated by a 
big bang at a later time. We call this the effective big bang. Of course, around 
this time we expect reheating effects to start to dominate, so that in reality 
the universe must pass over to a radiation-dominated era. However, the naive 
'effective big bang' concept is useful for extracting some qualitative predictions 
from our model. 

As well as giving a new effective start time for the big bang, inflation must 
also resolve the horizon problem. The conformal picture provides a straight- 
forward means of visualising this process. The key is to compute how much 
conformal time has elapsed by the end of the inflationary era. Using the same 
parameters as in figure |H1 we find that the conformal time evolves as shown in 
figure El As the universe emerges from the (real) big bang 77 initially grows at 
£ 2 / 3 , as can be seen from equation (|37|l . This gives rise to the rapid growth seen 
in figure which plots 77 as a function of ln(t/i p ). But once the inflationary 
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Figure 8: Hubble function exit from the inflationary regime. As the inflationary 
era ends the density and pressure start oscillating around values for a matter- 
dominated cosmology. The effective singularity for this dust cosmology is dis- 
placed from t = Q. The broken line plots the natural log of 2/3(i — 14.56). The 
parameters in this model are bo = 2.48, 64 = —0.51 and m = m p . 
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Figure 9: Evolution of the conformal time rj(t) as a function of ln(i). The 
parameters for the model are as given in figure |H1 and again t is measured in 
units of the Planck time. This plot sets fj, = 1, so the time variable must be 
scaled to correspond to more physical values of/iw 1CT 6 . 
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Figure 10: The horizon problem for a closed universe in the conformal picture. 
The lower solid line represents the surface of last scattering, and the broken line 
just above it is the current epoch. The thick line is the future boundary (t = oo). 
Null geodesies lie at 45°, so there is no problem with points in the observable 
part of the surface of last scattering having once been in causal contact. 

region is entered, R(t) starts to increase rapidly. So the conformal time, which 
involves the time integral of R , quickly saturates, as can also be seen in fig- 
ure[5] Any further increase in n will take an extremely long time to occur, since 
the integrand 1/R is extremely small. So if n has not reached tt/2 before R has 
inflated significantly, the universe will have to exist for an extremely long time 
to reach the boundary value of tt/2. 

For the parameters used in figures |H1 and El we can see that n saturates at a 
value of around 0.97 by the end of inflation (at around t = 100i p ). The value 
of R at this time is 2 x 10 22 /i _1 Z p . All changes in n beyond the end of inflation 
occur very slowly. Irrespective of the properties of the reheating phase, n will 
have changed very little by the time recombination occurs. This pushes the 
era of recombination to a much later conformal time, close to the value for the 
present epoch. This removes the horizon problem, as the obervable patch of the 
universe has been in causal contact and has had sufficient time to thermalise 
(see figure irn|) . This is not true of the universe as a whole, as antipodal points 
remain outside each other's light cone until t = oo. 

In order to make some initial estimates from our boundary condition, we 
will ignore the effects of reheating and nucleosynthesis, and assume that the 
universe just smoothly enters a matter-dominated phase described by 



For the parameters we are currently considering we have R% = 2 X 10 22 / p , 
to = 14.56i p and t\ = 100i p . We will also ignore the cosmological constant, 
though this will clearly affect the behaviour at late times. The conformal time 
elapsed at some later cosmic time t can be approximated by 




(41) 



/3 (h - t Q ) 2/3 « 0.97 + 2.4 x 10~ 19 i 1/3 , (42) 



l 



16 



0.2 0.4 0.6 0.8 1 

Figure 11: Matter-dominated trajectory in the J7m~^a plane. The trajectory 
(solid line) is picked out by demanding that the total conformal time elapsed 
throughout the evolution is 0.6. This evolution joins onto an inflationary model 
which takes up around 0.97 units of conformal time, so that the boundary 
condition of a total of tt/2 is obeyed. 



where t is measured in units of the Planck time. We can therefore arrive at a 
crude estimate for the age of the universe simply by setting rj — tt/2. This gives 
a time of 1.6 x 10 61 in Planck units, or about 2.8 x 10 10 years. This is a good 
illustration of how some very large numbers (such as the age of the universe in 
Planck times) can emerge from a simple model of inflation, augmented with our 
closed universe boundary condition. 

In order to produce a more detailed prediction, we return to considering 
trajectories in the Qm~&A plane. We can still assume that the universe is 
essentially matter dominated throughout its history, but now the total amount 
of conformal time that should elapse along the trajectory is less than ir/2. The 
effect of this is to drive us onto trajectories which are closer to the critical 
line. If we only require a total conformal time of around 0.6 to elapse during 
the matter-dominated phase of the universe (as for the model above), then we 
require an a value of around 10 4 , producing a trajectory that is very close to 
spatially flat, and predicting a universe which is currently closed at around the 
1% level (see figure ITT|) . 

We can now see that our boundary condition gives rise to a remarkable see- 
saw mechanism, with two key components. The first component relates to the 
amount of inflation that occurs. For a range of conditions, the time at which 
inflation starts can be approximated by 

A check on the physical approximations we make here is that this quantity has 
the correct scale transformation properties under the transformation of equa- 
tion (|36l) . The value of H at the start of slow-roll inflation can be approximated 
by l/(3ui), 

Hi = ^~, (44) 
Su, 
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and a suitable condition for the end of slow roll can be taken as 

,2 



H» f = fL. ,45, 



During slow roll the evolution of H is well approximated by a linear fit, so 

H = H i -^—^{H i -H f ). (46) 

Uf-Ui 

The gradient is determined from the field equations to be 

Uf — Ui 3 

The total number of e-foldings, N, can therefore be approximated by 

Numerical simulations reveal that this approximation does give a correct order- 
of- magnitude for the number of e-foldings, though it does tend to underestimate 
the precise value. As expected, the result for TV is invariant under scale changes, 
as is clear from the integral formula for N. The key conclusion is that the number 
of e-foldings is primarily determined by bo, with the value of H at the end of 
inflation fixed by p. The value of 64 essentially decouples from N. An increase 
in bo leads to an increase in the final value of R, and hence to an increase in pR 3 
in the matter phase. From equation l|2U|) we see that, in order to preserve a, 
any increase in pR 3 must be offset by a decrease in the cosmological constant. 
Conversely, the observed value of A can be employed to constrain the amount 
of inflation. 

The second quantity of interest, given our boundary condition, is the amount 
of conformal time that has elapsed before the universe enters the present matter- 
dominated epoch. Figure shows the typical evolution of rj. We see that a large 
percentage of the conformal time is taken up before inflation starts. During 
this region we can approximate R by equation 1)3 7|l. Taking Ui as roughly 
approximating the point where r\ starts to saturate, we find that the conformal 
time that has elapsed is approximated by 

And again, the result naturally assembles into a scale-invariant quantity. This 
result typically underestimates the total elapsed r\ at the end of inflation, but 
it does set a useful restriction on the order of magnitude of the expression on 
the right-hand side. We now see that, for fixed bo and p, the greater the initial 
curvature (increasing I&4I), the greater the amount of conformal time taken up 
during inflation. This in turn forces the universe today to be closer to flat. 
And conversely, the more spatially curved we want the universe to be today, 
the closer to flat it must be initially. This is the second component of the see- 
saw mechanism. The observed values of the cosmological parameters, together 
with our boundary condition, therefore place a severe restriction on the allowed 
models. 
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Given the strong constraints our model places on trajectories in the J1m-!1a 
plane, the model must predict the small value of the cosmological constant. To 
see how small values of A naturally arise, we return to equation l|20(l and write 

27a , . 

A -(8^W (50) 

In order to make some useful approximations we will again ignore the effects 
of reheating and radiation domination, and simply assume that p m R 3 is given 
by its value for the scalar field at the end of inflation. If we furthermore ignore 
the effects of curvature and A at this stage, we can write 8nGp m = 3H 2 , so we 
have 

3a 



A 



(51) 

end 



This expression does not depend on the precise moment of evaluation. As we 
have seen, typical values of a are in the range 10 3 to 10 5 . The small value of A is 
therefore seen to be a consequence of the large values of RH 2 ^ 3 achieved by the 
end of inflation. The value of H at the end of inflation is given by equation , 
and for R at the end of inflation, Rf, we can write 

R f = e N Ri. (52) 

Here Ri is the value of R at the start of inflation, and N is the number of 
e-foldings. We can estimate Ri using equation 1(37(1 and the value of Ui from 
equation 1)43(1 . Putting these together, we arrive at the approximate formula 

A « 8 x l0 4 aM^V 6J V. (53) 

The typical models we consider have N in the range of 40 to 50, and \i around 
10 -6 , which does indeed predict an extremely small value for the cosmological 
constant. As commented in the introduction, the dominant behaviour in A is 
governed by the exp(— 6iV) term, which immediately yields values for A of the 
expected magnitude. 



5 The curvature spectrum 

Our proposed 'quantisation' condition on the available conformal time places a 
strong restriction on the class of allowed models. We must now see whether any 
of these models produce a perturbation spectrum consistent with the observed 
primordial fluctuations in the CMB. Due to the scale-invariant nature of the 
equations (under the transformation of equation l)28[l). the parameters that 
determine the key features of the inflationary period are the scale-invariant 
quantities bo and &4//1 4 / 3 . Furthermore, the restriction that the total amount 
of conformal time available is it/2 constrains the combination in equation i|49|) 
to be of order 1. This limits us to a highly restricted class of models. 

The scale transformation we have utilised at various points is not conserved 
by quantisation, so vacuum fluctuations impose an absolute scale on the model. 
The key quantity that determines the size of the perturbations is the mass ratio 
fi = m/rrip. For models with a quadratic potential to generate perturbations 
consistent with observation we need to set ji to be around 10 -6 . As explained 
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above, this value is irrelevant for the evolution of the scalar field, but /x does 
set the absolute size of the Hubble function, and so fixes the magnitude of 
perturbations. 

In order to consider a concrete example, we adopt the following set of pa- 
rameters: 

6 = 2.37, 4tt = -0-49, M = 2.1 x 1CT 6 . (54) 

These values are typical for giving good agreement with experiment for the 
curvature spectrum. Evidently, little fine tuning is required in chosing the two 
scale-invariant combinations, both of which are of order unity. Furthermore, 
the predictions we make here are quite robust under variations of these pa- 
rameters. We are a long way removed from areas of parameter space where 
chaotic dynamics of the type considered by Cornish & Shellard ^B] could prove 
significant. 

In a closed universe model there are two scales of interest: the horizon size 
l/H, and the radius of the spatial sections R. The ratio of these defines the 
dimensionless quantity which is related to the closure density by 

°- 1= (W< (55) 

where the total SI is defined by 

n = n A + n M . (56) 

During inflation Q should decrease, as the universe is driven towards flatness. 
This effect can be seen for our chosen parameters in figure IT2l which confirms 
that the universe is accelerating over the range e 11 < t/t p < e 16 . The plot also 
shows that the radius is always greater than the horizon size, which is important 
when considering perturbation modes. 

The evolution of the scale factor for our chosen parameters is shown in 
figure H31 During the inflationary region ln(i?) increases from about 14 to 60. 
This corresponds to a total number of e- foldings of 46, as expected given our 
choice of parameters. This is roughly the figure we expect for the number 
of e-foldings if the universe departs from flatness by an amount that is visible 
today jS]- In order to give rise to the observed fluctuation spectrum it is thought 
that around 40-50 e-foldings are required between the point where the present 
horizon scale leaves the horizon and the end of inflation [T§|. For this to hold, 
we therefore require that the present horizon scale leaves the Hubble horizon 
soon after the onset of inflation. 

Closed universe models have an advantage over flat models when discussing 
scales and perturbations, because equation l|55|l gives a direct expression for 
R in terms of f2 and H. To apply this, suppose that at the current epoch 
a given physical size occupies a fraction x of the Hubble horizon 1/Hq. The 
corresponding physical scale is 

d = < x < 1. (57) 

no 

For wave modes, do will be quantised in units of 2TrRo/(n 2 — l) 1 ^ 2 , where n is 
an integer. We will ignore this quantisation effect here, as the number of modes 
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Figure 12: Evolution of 1/(RH). The natural log of the dimensionless ratio of 
the Hubble radius to the curvature radius is plotted against log time (in units 
of t p ). The parameters for the model are as described in the text. During the 
inflationary period, \/{RH) is a decreasing function of time, corresponding to 
an accelerating phase. 



60- 
50- 
40- 

ir 

30- 
20- 




10 4, - — ; ~ ~ , , , , 

6 8 10 12 14 16 

ln(t) 

Figure 13: Evolution of the scale factor. The evolution of logarithm of the scale 
factor R/lp is shown as a function of log time (t in units of t p ). 
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Figure 14: Evolution of the total CI during inflation. The total density is quickly 
driven close to flat, after which CI only changes slowly. 



available ensure that a continuum is a reasonable approximation for all but the 
low-fc modes. Now consider the distance corresponding to do m a- n epoch when 
the scale factor is R instead of Rq. This scale is simply Rx/{RqHq). This 
distance is equal to the Hubble horizon at the epoch corresponding to R when 

RX 1 (58) 



RqHq h 

(As is standard practice we formulate the conditions for horizon crossing in 
terms of the Hubble horizon, as opposed to the particle horizon.) It follows that 
the scale corresponding to a given x crosses the Hubble horizon when 



x = Ml = ^El (59 ) 

So as CI is driven to unity during inflation, progressively smaller scales cross the 
horizon as time advances. 

Figure ITT1 shows the behaviour of CI in our scalar field model. As expected, 
CI is driven close to 1 in the region e 11 < t/t p < e 16 . In order to use this plot in 
conjunction with equation (|59|l we need to assign a value to Clo. The most recent 
results from WMAP, with all data folded in, give a value of Cl = 1.02±0.02 [Tj, 
so here we work with a value of CIq = 1.02. The scale corresponding to the 
present horizon size (x = 1) is therefore obtained during inflation when CI = 1.02. 
Figure ITU shows that this occurs at \n(t/t p ) — 12.0. By reference to figure l*Hfl 
we see that this time corresponds to roughly 2 e-foldings. This leaves roughly 44 
c-foldings to the end of inflation, which are sufficient to produce a fluctuation 
spectrum with the desired properties. These values are consistent with the 
estimates of Uzan, Kirchner & Ellis. 0. 
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Each value of x less than 1 corresponds to a scale which leaves the Hubble 
horizon at successively smaller values of f2. A remarkable feature of the approach 
outlined here is that it does not require any knowledge of the evolution of the 
universe to relate scales today to scales at the end of inflation. This is because 
in any closed (or non-flat) universe it is possible to measure R directly. This 
cannot be done if the universe is assumed to be flat throughout inflation, as one 
loses the ability to track scales through the evolution of f£. In flat models the 
only way to find the scale ratio R/Rq is to evolve the scale factor R through 
the history of the universe, including the epoch of reheating, about which little 
is known. Of course, our method would be impossible to apply if the value of 
Slo turns out to be undetectably close to 1. This is predicted by approaches 
such as chaotic inflation, which can typically allow of order 10 10 e- foldings |19| . 
But for models that allow FIq to deviate measurably from 1 the method wc 
have outlined here should be of considerable use in calculating the fluctuation 
spectrum induced by inflation. 

In order to compute the primordial curvature spectrum we adopt the slow- 
roll approximation, following the presentation of Liddle & Lyth \F3\. For a 
quadratic potential, the slow roll condition can be expressed simply in terms of 
the single parameter e, where 

' (60) 



4lT(f> 



For slow roll to hold we require that e< 1. However, the model we consider here 
has e in the range from to 10 _1 , with the larger values corresponding to 
higher k values. We are therefore close to the point where slow roll ceases to be 
a good approximation. A concrete example of this is seen when we consider the 
spectral index (see also Martin & Schwarz (20j for a discussion of the precision 
of slow roll). For the main evolution of <fi, however, plots such as figured 
confirm that the slow roll approximations are holding up, with the evolution 
of H roughly linear with the expected gradient of /x 2 /3. A further source of 
error, as stated earlier, is that we are ignoring the quantisation effects implied 
by a closed cosmology. But, despite these provisos, we may hope that the 
simple, slow-roll equations will still give a reasonably accurate prediction for 
the curvature spectrum. 

Given the slow roll approximations, the key equation determining the curva- 
ture spectrum Vn(k) is 

™>-(f)'(£)"' < 6i » 

where the right-hand side is evaluated when a given scale x crosses the horizon. 
If we write 

= (62) 

27T0 

we see that the power spectrum is controlled by the ratio H 2 /<f>. This function 
is plotted in figure ll5l Figures ITU and lTBl contain all of the information required 
to extract V-jiik). Given a comoving wavenumber k, we find x from 

k H a 
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Figure 15: Evolution of H 2 /(f) during inflation. The function plotted is 
— 10 5 H 2 /cf> as a function of \n(t/t p ). This function controls the magnitude of 
the curvature perturbation. 

The value of x is converted to a time using equation (|59|l and the relation- 
ship plotted in figure El With the value of t then found, we can read off the 
spectrum from figure El The fact that the graph turns over for ln(t/t p ) < 12 
suggests that the spectrum will contain less power at low-fc values than would 
be expected from a straightforward power law in a flat cosmology, in contrast 
to the predictions of Starobinsky 5 . But again, more detailed calculations are 
required to confirm this point. 

It is convenient to plot perturbation spectra as a function of y = ln(l/a;). 
With dimensions reinserted, y is given by 

y = In ^ x 3 x 10 3 ^ (64) 

where fc is measured in Mpc -1 and the Hubble constant is given by Hq = 
100/ikms -1 Mpc -1 . In figure [TBI we plot the power spectrum as a function of y 
for our chosen parameters, including flo = 1-02. (For a plot showing the power 
spectrum as an explicit function of fc in Mpc -1 see figure El) The quantity 
plotted (solid curve) is 10 lc '(27r) 2, P-R,(fc) ■ The graph shows a clear cutoff at low- 
fc values. A cutoff of this type was proposed on phenomenological grounds by 
Efstathiou [5JJ • Indeed, it turns out that the type of exponential cutoff proposed 
by Efstathiou agrees remarkably well with our calculations, as can be seen in 
the dashed line, which plots 

10 5 2nT 1/2 {k) = 37.6(1 - .023y)(l - exp(-(y + 0.93)/0.47)), (65) 

valid for y > —0.93. The fact that our model predicts a spectrum which has 
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Figure 16: The power spectrum of curvature perturbations. The function plotted 
is 10 1C \2n) 2 VTi{k) as a function of y = ln(l/x) = ln(fc) + constant. The solid 
line represents the numerical predictions from our model, and the broken line is 
the best-fit power law with an exponential cutoff. The parameters are described 
in the text. 
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already been argued for on phcnomenological grounds is very reassuring. We 
can now see that our model will produce a deficit in the CMB power spectrum 
at low t, and in section [S] we confirm that this is the case. 

FigurefTBlcontains a further surprise. From equation l(6"5|l we see that at large 

1/2 

k the relationship between (k) and ln(fc) is roughly linear, which is quite 
different from the expected power law relation. To understand why this is the 
case we return to the slow roll equations and derive the form of the spectrum for 
a quadratic potential, following section 7.5 of Liddle and Lyth [T§]. The main 
relation we require relates the derivatives with respect to k and <fi: 

" ' " (66) 



dink Airtfi dtp 

where we specialise to results for a quadratic potential. During slow roll the 
following approximate relations hold: 

H 2 = ^ ftp, h=^=- (67) 



It follows that the curvature spectrum is approximated by 



We therefore find that 



K (k) = 4(-) [i<i> ■ (68) 



^V]i\k) = -%. (69) 



dink n v^F' 

1 /2 

This confirms that the relationship between (k) and ln(fc) is indeed linear. 
Furthermore, the gradient is predicted to be — 2///V37T. For the model we 
are considering this evaluates to 1.35 x 10~ 6 , which agrees well with the fit of 
equation (gSJl, which has a gradient of 37.6 x 0.023/(2tt x 10 5 ) = 1.37 x 1(T 6 . 
We can also see that 

H V lz (k)) = -2-±- ^ = -_L = _4e (70) 



dink p^ 2 (k)V^ W> 2 

where e is the slow-roll parameter of equation (|60|) . The typical slow-roll ap- 
proximation is to now write V-jiik) oc k~ , so that 

V^ 2 (k) fts Ak- 2e k A(l - 2eln(fe) + 2e 2 ln 2 (fc) + •••)• (71) 

The higher order terms here are incorrect, and introduce an error of order e ln(k) 
to the power spectrum. As commented earlier, at high-fc values e is in the range 
0.01-0.1, so this deviation can become quite significant. 

In figure IT7I we compare our computed scalar power spectrum with the best 
fit power law (n s = 0.96) and the WMAP running spectral index best fit. The 
vertical normalisations for these fits have been chosen arbitrarily, since it is the 
shape of the spectrum at large k which is of greatest interest here. The WMAP 
running spectral index model shape is given by 

V n (k) = Vn(ko) exp ({n s - 1) bx(k/k ) + \n rm (ln{k / k )f\ . (72) 
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Figure 17: Comparison of the scalar power spectrum with power law models. The 
function plotted is 1Q 10 (2n) 2 Vii(k) as a function of log 10 (fc), assuming h = 0.65. 
The solid line represents the numerical predictions from our model. The long 
dashes represent the best fit power law (n s — 0.96) and the short dashes are the 
WMAP running spectral index best fit. Notice that the vertical scale runs from 
800 to 1300, so the differences are slightly exaggerated. 



Using the WMAP figures for combined data sets their best fit values are n s = 
0.93, n run = —0.031 and ko = 0.05 Mpc -1 £Q. The running spectral index graph 
is interesting in that it suggests that this model is attempting to emulate both 
the cutoff at low k, as well as the reduced power at large k. The latter is possibly 
required as the strict power law spectrum at large k appears to over-predict the 
abundance of dwarf galaxies. This is discussed further in sectional 

A further quantity of interest is the tensor spectrum, which is significant in 
searches for B-mode polarisation in the CMB. As our definition of the tensor 
spectrum we take 

16 

V grav (k) = — H 2 , (73) 

IT 

where again the right-hand side is evaluated as the scale controlled by x crosses 
the horizon. (Here we follow the convention of Martin & Schwarz [201 )■ The 
main quantity of interest is the ratio of the tensor and scalar modes r, 

where the right-hand side is evaluated at some suitable low k. Following the 
slow-roll approximation, we can write 

, s 64uV , 
TWfc) = "IT- (75) 
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Figure 18: Power spectrum of tensor perturbations. The quantity plotted is 
10 10 (27r) 2 7 : 'giav(fc) versus y = ln(fc) + constant. The solid line is the numerical 
predictions. The dashed line is the best fit power law with exponential cutoff, 
offset by 20 units vertically for clarity. All parameters are described in the text. 

Applying d/(d Ink) to this now yields 

dnVk V — {k) = -3T' (?6) 

which is a constant. For the tensor mode spectrum we therefore expect that 
■Pgrav should be linear in ln(fc). 

Figure 1181 shows the predicted tensor mode perturbation spectrum for our 
model (again with fin, = 1-02). Also shown (slightly offset) is the best fit model 
of the form of equation (|65[) . together with an exponential cutoff at low k. This 
best fit is given by 

10 10 (2tt) 2 "Pg r av (fc) = 255(1- .023 y)(l-exp(-(y + 2.6)/0.65)) (77) 

valid again for y > —0.93. The slope in this formula is 255 x 0.023/(27r x 
10 5 ) 2 = 1.48 xl0~ n , which agrees well with the theoretical value of 32/Lt 2 /(37r) = 
1.47 x 10 -11 . If we form the ratio r at y — 2, which corresponds to k = 
0.0016 Mpc -1 for Hq = 65 kms -1 Mpc -1 , we find r ~ 0.19. This is large enough 
to be detectable in the near future. For example, the predicted sensitivity 
level to B modes for the Planck satellite corresponds to a detection limit of 
r > 0.05 |22|. comfortably below the value predicted here. A relatively large 
tensor mode component is therefore a firm prediction of our model. 

We end this section with a discussion of how to overcome the limitations of 
the slow-roll approximation and compute the primoridal power spectrum more 
accurately from first principles. The starting point for calculating the scalar 
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curvature spectrum in a simple flat model is to write the linearised perturbed 
action in the form |2~3*] 



S 2S= - 




x I v' 2 — rf-'v t iV i j H v 2 



) 



(78) 



where v is a gauge-invariant combination of matter and metric perturbations, 
and z is defined by z — <t>' /H. Here 4>q is the unperturbed field, and the dashes 
denote derivative with respect to conformal time. By working in this way the 
entire problem is reduced to one of analysing a scalar field with a time-dependent 
mass term. 

This approach generalises to the non-flat case, although here the definition of 
z becomes more complicated [2]. Also, for closed models, the mode expansions 
necessary to carry out quantisation have to be performed using spatial sections 
given by the 3-sphere S 3 . In this case the 'comoving wavenumber' k takes on 
integer values, with k = 3 its lowest (non-gauge) value. The mode equations 
found this way are very complicated, but retain the general form 



In the flat case f(rj, k) would be z" /z, and in all cases f(rj, k) is calculable from 
knowledge of the background evolution. 

The challenge now is to find suitable 'quantum initial conditions' so that after 
evolution through inflation, the variables Vk can be used to find the perturbation 
spectrum. To achieve this involves a mode decomposition of Vk into positive and 
negative frequency modes. The standard way to approach this is to assume that 
in the asymptotic past the background is either Minkowski or de Sitter, so that 
one knows the correct vacuum to chose. But this is clearly inappropriate here, 
as looking back in time we encounter the initial singularity. 

The question of how to proceed in the absence of any asymptotic notion of 
the vacuum state has been widely discussed ^3 EH1 • One simple approach 
is provided by Hamiltonian diagonalisation, where on each time slice modes are 
decomposed into positive and negative energy states by the Hamiltonian oper- 
ator. But this technique tends to overestimate the particle production rate. A 
clear way to proceed was developed by Parker and Fulling [23 ESJ and intro- 
duces the concept of an adiabatic vacuum. The application of this to the present 
case is complicated, but possible numerically, and initial results from this ap- 
proach support a low-fc cutoff, although with effects less pronounced than we 
have found using the standard approximate techniques. This will be discussed 
in future work. 

6 CMB power spectrum and large scale 
structure 

We now turn to using our primordial power spectra to compute the measured 
fluctuations in the CMB. For this purpose we use the CAMB package \%7\, which 
takes as its input the curvature spectrum and various cosmological parameters, 
and returns the Cg values of the CMB power spectrum. For the cosmological 
parameters we take ficdm^ 2 = 0.135, £l\>h 2 = 0.0224 and the optical depth to 
reionisation r = 0.17. Together with SIq = 1-02 these correspond to the WMAP 



v^ + (k 2 -f( Vl k))v k = Q. 



(79) 
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Figure 19: Predicted CMB power spectrum. These predictions are for a model 
with Qq — 1.02, as discussed in the text. The experimental points shown with 
la error bars are the WMAP results. 

best fit values pQ. Additionally, we use a Hubble constant of 65kms _1 Mpc -1 . 
This is lower than the best fit value quoted WMAP of 72 ± 5 km s" 1 Mpc" 1 . 
This latter value is obtained for an assumed flat model. If one relaxes the 
assumption of flatness, then the CMB data is quite degenerate with respect to 
the values of Hq and FIq (see section 6.2 of Spergel et al. pQ for a discussion of this 
point). Our value of Hq = 65kms _1 Mpc -1 is compatible at la with virtually 
all current determinations of Hq. For example, the HST Key Project value 
is 72 ± 8kms _1 Mpc _1 j2E], whereas combined Sunyeav-Zel'dovich and X-ray 
flux measurements consistently favour a slightly lower value [2H] ■ A lower value 
of Hq was also suggested for closed models by Efstathiou [2U, an d our chosen 
value provides a reasonable fit with the WMAP power spectrum for Hq = 1.02. 
It is worth pointing out in this context that the WMAP simultaneous best fit 
values of Ho and flo actually provide a rather poor fit to their measured CMB 
spectrum! 

In fixing £Im and Hq we in turn imply a value of Together these pick out 
a unique path in the Ojf— Q,a plane, for which we can calculate a total elapsed 
conformal time. In order to satisfy our boundary condition we need to achieve 
a value of conformal time by the end of inflation of 0.975. The value of 64//U 4 / 3 
introduced in equation 1)54(1 was chosen to achieve this. With all the parameters 
now chosen, the computed CMB power spectrum is shown in figure The fit 
is reasonably good, particularly at low I where our model does show a fall off 
in power. 

A further test of our model is the predicted linear matter power spectrum. 
This is shown in figure EH where it is compared with the 2dF data of Percival 
et al. .30;. Also plotted is a fiducial model with n s — 1, ag — 1, h = 0.7, 
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Matter Power Spectra 




o 



0.01 0.1 1 

k (h Mpc)~ 1 

Figure 20: The linear matter power spectrum. The solid line is the prediction 
from CAMB using as input our scalar primordial power spectrum. Long dashes 
are the 2dF measurements [30] convolved with the survey window function. 
The dot-dash line is a comparison fiducial model with n s = 1, erg = 1, h = 0.7, 
£l m h = 0.2, and = 0. 
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il m h = 0.2, and fif, = (and therefore no baryon oscillations). This model is 
calculated from the fitting formulae of Eisenstein & Hu ■ In order to achieve 
a good fit with the data, the vertical scale for our model predictions has been 
adjusted downwards by about 19%. Clearly, the model we have used so far needs 
some adjustment in its chosen parameters to correct the normalisation. But here 
we only we wish to compare the overall shapes of the three matter power spectra, 
leaving the overall normalisation aside. It is clear that the spectrum based on 
our non-power law primordial spectrum fits the 2dF data with a similar degree 
of fidelity to the power law-based fiducial model. On large scales (small k) both 
appear to overestimate the 2dF results. This can be explained by the fact that 
the 2dF results are given convolved with the survey window function, which has 
its largest effects at large scales (sec figure 2 of Pcrcival et al.). On smaller scales 
our results fit the shape of the 2dF spectrum very well. There is therefore no 
apparent problem in having a non-power law primordial spectrum right down 
to these scales. On even smaller scales the relevant comparison would be with 
Lyman-a data, but this area is still controversial and will not be discussed here. 

So far we have developed a model based on the best fit WMAP values for 
fi c d m /i 2 , fl^h 2 and Slo = 1-02, but using a different value of Hq. This model 
does provide a reasonable fit, but figure^] clearly leaves room for improvement, 
as does the normalisation of the matter spectrum. Ideally we would employ an 
MCMC analysis to find the best fit parameters for our model, but here we simply 
give an alternative model which improves the fit. This model has fio = 1-04, 
so it still consistent with WMAP at the la level. We also take fl b h 2 = 0.0224 
(as before), h — 60 and fl c a m h 2 = 0.110. Together these yield a value of 
fi c d m = 0.306, which is reasonable. The values of bo, 64 and /1 consistent with 
these are 

6o = 2.3310, -4k = -0.2798, M =1.94xl0~ 6 . (80) 

The lower value of \b±\ produces a universe with larger curvature today. The 
results of this model are summarised in a best-fit curvature power spectrum 
given by 

10 5 2nV^ 2 (k) = 35.02(1 - .023 3/) (1 - exp(-(y + 0.20)/0.397)). (81) 

The result for the CMB power spectrum with these parameters is shown in 
figures PHI and 1221 The former includes the WMAP points, and the latter plots 
both WMAP and the latest Very Small Array (VSA) data points, plotted on a 
linear scale • This model clearly gives an extremely good fit to the observed 
spectrum, with the dip at low £ now quite pronounced. This is because an 
increase in f2o has the effect of sliding the graph to the right in k space, as is 
clear from equation l|59|l . There is also a good fit with the VSA data at higher 
I values. The main disadvantage with this model is the low value of Ho, though 
this is still within 1.5a of the HST Key Project value. 

7 Discussion 

Many authors have argued that it is difficult to obtain realistic models of a 
closed inflationary universe 0J El El Ej ■ One of the main problems is the 
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Figure 21: CMB power spectrum for a model with Qq = 1.04 (I). The parameters 
are discussed in the text. The experimental points shown with la error bars are 
the WMAP results and the dashed curve corresponds to the best fit A CDM 
power law CMB power spectrum as distributed in the WMAP data products. 

Cg predictions 
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Figure 22: CMB power spectrum for a model with Qq = 1.04 (II). A similar plot 
to figurc[^ but now plotted on a linear scale, with the latest Very Small Array 
data points also included to extend out to higher I values (open circles) . 
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low number of e-foldings required, which it is argued demand an unacceptable 
degree of fine tuning in the initial conditions for inflation. We have shown that 
this is not the case. It is relatively straightforward to construct inflationary 
models resulting in a universe which is measurably closed today, using only the 
simplest model of a free massive scalar field. These models have a number of 
attractive features and provide a very good fit to the observational data. 

Motivated by the global geometry of de Sitter space, our model incorpo- 
rates a quantisation condition on the available conformal time. This may seem 
uncomfortable from the point of view of causality, but it is quite sound from 
the viewpoint of wave mechanics. In field theory one regularly applies bound- 
ary conditions at temporal infinity, and the model presented here extends this 
idea to the universe. Furthermore, the constraint does not relate to the actual 
conformal time elapsed in the universe, but to the available conformal time at 
the point of creation of a homogeneous, idealised universe, before substantial 
structures have formed. 

We hope that this quantisation condition will serve as a guide in constructing 
more detailed physical models of the processes taking place around the big bang, 
and that it might emerge from a quantum gravity understanding of this phase. 
Of course, many quantum gravity models of the very early universe involve 
tunneling from a Euclidean phase, in which case one cannot argue against the 
constraint on grounds of causality. Furthermore, for those who remain unhappy 
with our proposed constraint, we would point out that many of our conclusions 
do not depend crucially on it. One can relax the condition, explore the full 
space determined by the input parameters, and still make sensible predictions 
about the primordial power spectrum. 

With a massive scalar field included, we find that initial data can be specified 
in terms of a pair of scale-invariant coefficients in the expansion of the scalar 
field around the initial singularity. These then determine the initial conditions 
for inflation. Quantum gravity effects will not be significant until well before the 
onset of inflation, so we are fully justified in running the evolution equations 
back to much closer to the initial singularity. Given only that we require of 
the order of 44 e-foldings we can conclude that the scale-invariant parameters 
bo and 64//J, 4 / 3 are both of order unity. Because we start close to the initial 
singularity, with infinite if, the model predicts a period of inflation in the early 
universe, independent of any question of generating super-horizon fluctuations. 
Furthermore, our condition on the total conformal time ensures that the matter- 
dominated plot in figure|H|represents the largest departure from flat that we can 
consider. The inflationary period drives us closer to flat and, assuming f2jvf ~ 0.3 
we are led to an Da of around 0.7. Viewed in reverse, given input parameters 60 
and &4//1 4 / 3 of order unity, together with our boundary condition, we predict an 
extremely small value for A, consistent with observation. The final parameter 
in the model is the mass ratio fi = m/m p , which is fixed by considering the 
growth of perturbations during the inflationary period. 

The predictions of our model are robust against small variations of the input 
parameters. There is no question of the initial data being specified in a chaotic 
region of the classical dynamics. Chaotic behaviour in the inflationary equations 
in a closed universe was noticed by Page and Cornish & Shellard ^Hj- But 
these authors considered initial data around H — 0, whereas our parameters 
specify the fields as they exit the initial singularity, with H = 00. The evolution 
therefore never enters the region of parameter space where chaotic evolution 
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becomes relevant. The main element of fine tuning required in our model is 
setting ji = m/m p to I0 -6 , but this type of tuning of the potential is a well- 
known feature of all inflationary theories ^7]. Masses of the order of 10~ 6 m p 
are quite legitimate from the viewpoint of high-energy physics, which we expect 
will shed light on the possible origins of such particles. 

Our model makes a series of predictions. The universe is closed at around the 
level of a few percent, and we can explain the dip in the CMB power spectrum 
observed at low I. The statistical significance of this dip has been discussed 
by many authors [HI El ESI EE! • Their 

main conclusion is that this dip is not 
inconsistent with the spatially flat concordance model, but could still be an 
indicator of new physics, provided that any new model does not introduce too 
many additional parameters. In this respect our model is very efficient. Taking 
Qjif and Q\ as input data introduces two constraints on the three parameters 
bo, 64 and ji. The one remaining degree of freedom is removed by fixing the 
scale of the overall perturbation amplitude. We also predict a running spectral 
index with reduced power at high-^ values. This looks to be favoured by dwarf 
galaxy and Lyman a data, though this area is controversial. The behaviour of 
the CMB spectrum beyond the second peak will also help break degeneracies 
in various models, and high-resolution measurements in this region would be 
extremely valuable. 

For our best-fit models we typically predict of the order of 50 e-foldings. Such 
values are quite natural in our scheme, removing the central objection to closed 
universe cosmologies. Furthermore, a value of 50 has the additional feature of 
removing any potential problems with trans-Planckian physics |17j . That is, 
there is no question of Planck scale physics being inflated up to cosmological 
scales. A further indicator that our predictions are relatively insensitive to 
quantum corrections is provided by equation Ij40|l . which gives the value of the 
total energy in the universe as we leave the big bang. Assuming that 64 //1 4 / 3 is 
of order unity, this expression becomes 



f ~ °- 03 h rs2i 



This tell us that the action Et is very large, when measured in units of h. This 
is reassuring, as it justifies treating the universe on the whole as a classical 
object. While quantum gravity effects will certainly alter how the universe 
behaves around the singular region, there is good reason to believe that the 
predictions of this paper will be largely unaffected by any such new physics. 
Furthermore, the limiting value of the action as t h-» is large in units of h, 
suggesting that the universe started in a highly improbable state. This will 
naturally favour spherically-symmetric initial conditions, and could also prove 
attractive from the viewpoint of the extremely low entropy that the universe 
has to have initially. 

One further aspect of our model where considerations involving energy are 
relevant is the end of inflation, as the universe enters a reheating phase. Con- 
verting the energy density at the end of inflation into a temperature yields the 
reasonable value of 10 16 GeV. This is as expected for typical reheating, if it 
occurs immediately after inflation is over. But as the scalar field exits inflation 
the universe is in a matter-dominated phase. Some process has to intervene 
which converts the scalar field into a radiation-dominated phase. Curiously, 
this may be difficult to achieve in a way which preserves the total, integrated 
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energy, though there is little difficulty in demanding local conservation. This 
is probably a reminder that the 'total energy' in a closed universe remains a 
problematic concept in general relativity. 

An important question is to identify the main experimental observations 
that will argue for or against our model. Our first key result is the low-^ dip 
in the power spectrum. The significance of this could be further enhanced by 
any lowering of the error bars in these measurements. While cosmic variance 
is the ultimate determining factor there is still some room for improvement, 
particularly in the question of foreground contamination M 1 . The drop in 
the power spectrum at large I is also a characteristic feature, which could be 
detectable with future observations. The third key prediction is the high value 

of i the ratio of tensor to scalar perturbations. This should also be observable 

by experiments such as the Planck satellite. It would also be of great interest 
to see which of these predictions is dependent on our choice of potential, and 
which survive if this restriction is removed. 

We hope that much of the future work suggested by this paper is self evi- 
dent. Before significant future work can begin we need to address the two main 
provisos mentioned in the text. The first is that mode quantisation in a closed 
universe has to be built in. The second is that the mode equations need to 
be integrated numerically, as aspects of the slow roll approximation are inade- 
quate for our model. Neither of these pieces of work should prove too difficult, 
and a start has been made by Zhang and Sun [21] who appear to support our 
use of equation l|tjl|l as a leading order approximation. One area that we aim 
to address quickly is to explore in more detail the effects of varying the main 
input parameters, to explore the range of possible universes expected today. 
One significant point here is that all evolution is constrained to lie inside the 
matter curve of figure El This in itself is a significant restriction. In addition, 
the parameters chosen for the examples presented here were not computed by 
detailed model fitting. A full Bayesian determination of these parameters using 
all available data is clearly desirable, and will give a more stringent test of how 
well our model holds up. 
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A Conformal line elements 

Here we establish spacetime conformal forms for the standard cosmological line 
elements for open, closed and flat cosmologies. (The fact that the Weyl tensor 
vanishes for cosmological models ensures that this is possible.) Our starting 
point is the FRW line element in the form 

ds* = dt* - (i + ^ 2 2 /4)2 (dr- + r\d^ + sin 2 W d^)). (83) 

This line element is already conformal in its spatial component. But a problem 
with this form is that r is assumed to be dimensionless. To rectify this we 
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introduce a constant distance A and replace the line element with 

ds " = df2 - (tf+kr'y ( dr2 + r 2 (de 2 +s\n 2 {6)d4> 2 )). (84) 

Here t, r, A and R all have units of distance (assuming c — 1). As usual, the 
constant k is either ±1 or zero. 

We seek a coordinate transformation to the spacetime-conformal line element 

ds 2 = (dt 2 - dr 2 - x 2 [d9 2 + sin 2 (6) d(f> 2 )) , (85) 

where / is a (dimcnsionlcss) function oft and r. Clearly, from the angular terms, 
we must have 

r _ 2XrR(t) 

/-^W (86) 

In addition, the coordinate transformation must satisfy 

r 

dt = f cosh(w) dt H — sinh(u) dr 
r 
r 

dr = f sinh(u) dt H — cosh(u) dr. (87) 
r 

For flat cosmologies (k — 0) we simply set t/A equal to the conformal time rj, 
with r = r and the hyperbolic angle u set to zero. For non-flat cosmologies it 
is perhaps surprising to find that u is non-zero. That is, there is a mismatch 
between the conformal coordinate frame and the cosmological frame. 

The integrability conditions for the coordinate transformation provide a pair 
of differential equations for u: 

Ou 

r— = sinh(u) (88) 
or 

2\rR(t)^ = (A 2 + kr 2 ) cosh(w) + (-A 2 + kr 2 ) (89) 

Equation l|88(l can be solved straightforwardly to give 

1 + rA{t) 

1 - rA(t) ■ 1 j 

The function A(t) is determined by equation (|89J) . which reduces to 

2XRA = \ 2 A 2 + k, (91) 

where the overdot denotes the derivative with respect to cosmic time t. The 
solution for A depends on the curvature: 

XA=H^ k = 1 (92) 
l-tanh(ry/2) k = -1. 

The constant of integration has been chosen such that A = corresponds to 
r/ = t = 0. 

Next we solve the following pair of equations for r: 

dr dr 

r— = rcosh(w), — = /sinh(u). (93) 
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With a suitable choice of the arbitrary scale factor we find that 

\ 2 A 2 /k + l 



1 - r 2 A 2 

The remaining equations determining t arc 



r. (94) 



dt dt 

r— = rsinh(u), — = /cosh(u). (95) 

These are solved by 

A 2 + kr 2 A 

t= i^T (96) 

Here the arbitrary constant of integration has been fixed to ensure that t = T] = 
corresponds to t = 0. 

With r and t now found, we can write 

/ = ^^" t (97) 

Our remaining task is to find an expression for t in terms of t and r. For this 
we use that result that 

2At _ 2XA 
X 2 + k(r 2 - t 2 ) ~ k - X 2 A 2 ' ^ ' 

The right-hand side is now a function of t only, so the left-hand side can be used 
in place of cosmic time. Again, the explicit form depends on the cosmology. For 
spatially closed cosmologies we can write 

2At 

A2 + r2 _ t 2 = tMv) (99) 

whereas for spatially open we have 

2At 



tanh(f7). (100) 



A 2 + t 2 — r 

The three cosmological scenarios can now be summarised as 

f = Rskn) = g U+^-t 2 ) I closed 

/ = ^ flat (101) 

t _/ 2At \ t 

/= i^mh^ = - 9 U 2 +t 2 -r 2 J A ° pen 

which arc the results used in the main text. 
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